The American Naturalist

Rapid evolution allows populations to persist in environments where they would otherwise go extinct. This phenomenon, known as evolutionary rescue, is typically studied in the framework of biological evolution, yet adaptive traits can also arise and spread through cultural evolution. The present study developed a stochastic eco-evolutionary model to compare rescue probabilities through biological and cultural evolution. Transmission bias governed the rescue probability under cultural evolution by setting how readily a rare adaptive trait was copied. Conformity bias suppressed population persistence because a rare trait was the least likely to be copied. Content bias toward the adaptive trait enabled evolutionary rescue when social learning was rapid, but it typically yielded a lower rescue probability than biological evolution. Only anticonformity bias, together with a high social learning rate, exceeded the rescue probability of biological evolution by enabling the adaptive trait to be established more rapidly. These results demonstrate that transmission bias alters the demographic consequences of cultural evolution and highlight the importance of transmission processes in evolutionary rescue theory. Understanding how adaptive behaviours are socially transmitted may also improve predictions of animal population persistence and inform conservation efforts in rapidly changing environments.

The shape of density-dependence governs species persistence, and ecosystem stability. Yet, whether per-capita growth declines sublinearily, or superlinearily with density remains hotly debated. Growth rates across the tree of life have been shown to decline sublinearly with density, whereas theory founded on resource competition predicts the opposite. Here, we resolve this discrepancy and show that sublinearity can readily emerge from geometric constraints on consumer interactions. By linking inter individual spacing, movement and interference rates, we derive two limiting-interference regimes, one of which the well-mixed limit recovers the form of classic Beddington DeAngelis interference response. We then developed an individual-based model from first principles which reproduces the derived sublinearity response, and further use empirical data from published consumer-resource experiments that also bears the signature of sublinear density-dependence. Further, embedding the interference mechanisms underlying the emergence of sublinear density-dependence in coexistence theory opens a new regime for species coexistence where classical theory fails to predict. Our framework indicates that non-consumptive interactions are not merely a correction to resource competition but might be a distinct axis along which diverse communities may potentially coexist.

Ecologists remotely track movement steps of animals (e.g., via global positioning systems) and use step selection functions to study the effect of environmental factors upon their movement decisions. Constructing such functions requires pairing each observed step with some number of unobserved but feasible comparison steps. Larger numbers of comparison steps generally yield better estimates but also incur potentially challenging computational demands. Thus, it is important to determine an appropriate number of comparison steps. No established guidance exists for this decision. Here, we use simulated tracks to assess how many comparison steps are needed, fitting each set of steps to a conditional logistic regression model. We monitor errors in estimated effects for several classes of tracks, identifying the number of comparison steps for which mean relative absolute error in estimated effects is consistently low. By this criterion, 32 comparison steps per observed step are needed for our primary class of simulated tracks. Tracks in more homogeneous landscapes, tracks with shorter mean step lengths, or shorter tracks generally require more comparison steps (ranging from 64 to 128 per observed step) to achieve the same level of accuracy. Longer tracks generally require fewer comparison steps (16 per observed step). These results clearly demonstrate that the number of comparison steps influences how well step selection functions estimate covariate effects and provides initial direction in a research area that currently lacks quantitative guidance. Movement ecologists should take care when selecting the number of comparison steps paired with each observed step because those decisions matter.

Heritable bacterial symbionts are pervasive in terrestrial arthropods, often imposing reproductive manipulations to promote their own spread within host populations. Co-infections are common, potentially allowing symbiont co-infectors to hitchhike through a host population. However, adverse thermal conditions can disrupt these communities, particularly when co-infectors vary in their thermal sensitivity. We used a multi-generation experiment to test whether warm (29 {degrees}C) conditions disrupted spread of heritable symbionts through uninfected populations of the spider, Mermessus fradeorum. We tested two common infection combinations: a single infection with a cytoplasmic incompatibility (CI) inducing Rickettsiella or a feminizing co-infection that included a feminizing Wolbachia, the same Rickettsiella, and up to three apparent hitchhikers (two additional Wolbachia strains and Tisiphia). We initiated replicate populations with 1/3 of one infection type and 2/3 uninfected spiders, evaluating population infection rate over 5 spider generations under different temperature regimes. Under cool (21{degrees}C) conditions, Wolbachia feminization drove co-infection to 88% and Rickettsiella CI drove single infection to 83% of host populations. Vertical transmission for all symbionts was high (97-99%) and hitchhiking symbionts also spread effectively. Under warm conditions, feminization and CI efficacy were reduced, and symbionts suffered variably reduced vertical transmission. Warm conditions ultimately destroyed the co-infecting symbiont consortium and impeded symbiont spread. On its own, though, Rickettsiella was still able to increase, despite reduced strength of CI. We hypothesize that contrasting tensions between feminizing spread of the symbiont consortium versus environmentally driven loss of function and transmission may explain observed patterns of mixed infections in field populations of this spider.

Squid abundance has been reported to increase globally between 1970 and 2010. This increase has been hypothesized to result from two primary factors: the loss of top predators due to overfishing and rising ocean temperatures. The decline in apex predators may lead to the expansion of squid populations either through reduced predation pressure or diminished competition with juvenile predators. Concurrently, increased temperatures could enhance the somatic growth rates of squid, thereby accelerating their population growth. However, empirically disentangling the impacts of predator loss and temperature on squid biomass remains challenging, especially in a food-web context. In this study, we used a size- and trait-based model of upper trophic levels that resolves the ecosystem structure -- biomass and trophic interactions of fish and squid -- for varying depth, temperature, and secondary production, to investigate two hypotheses of the historical expansion of squid, i.e., the effects of predator depletion from fishing and rising temperatures on squid biomass. Our model reveals that intensified fishing of squid predators -- specifically large demersal fish in shelf systems and large pelagic fish in open oceans -- leads to a slight increase in squid biomass. Conversely, elevated temperatures are associated with a decline in squid biomass. This temperature-driven reduction in biomass is attributed to an increased metabolism of squids beyond the available food supply. If historic overfishing on large marine predators continues to be curtailed, we expect a corresponding reduction in global squid biomass and fisheries potential, which could be further exacerbated by rising temperatures.

Across different modalities, signals play a core role in attracting mates and influencing mating success. In several non-human primate species, females produce calls during mating that are thought to promote male competition over receptive females. The extent to which social system characteristics modulate the function of copulation calls remains less clear. We studied copulation calls in wild Guinea baboons (Papio papio), who live in a multilevel society structured around units in which females associate and mate almost exclusively with a single male. We hypothesised that females use copulation calls as an indirect form of mate competition, with competition increasing in larger units. In addition, we hypothesised that females are more likely to mate again after calling. We analysed 6116 copulations between 2014 and 2025, involving 99 reproductively active females and 78 subadult and adult males. Females produced copulation calls in 72.7% of copulations, with large inter-individual variation. Neither unit size nor its interaction with the female's swelling size or the presence of simultaneously receptive females affected the probability of calling. A survival analysis with a subset of the data (2353 copulations) revealed no effect of calling on the latency to the next mating. Our results render the hypothesis that female Guinea baboons use calls in indirect mate competition unlikely. Yet, the probability of calling varied with sexual swelling size, suggesting that calls signal female fertility. Possibly, Guinea baboon copulation calls represent an evolutionary remnant, no longer under selective pressure, and can be considered index signals of female fertility.

1. Forest canopies harbour strong microclimatic gradients that shape plant performance, species interactions and ecosystem processes. Yet, despite renewed interest sparked by global change, forest canopies remain difficult-to-access experimental spaces. 2. With the goal to expand access to tree canopies as experimental arenas, we designed, built, and tested TreeTOP, a standardized experimental platform that opens canopy space for manipulative ecological experiments, specifically with potted plants. TreeTOP features lightweight aluminum frames placed in mature tree canopies non-invasively, allowing potted plants to be placed in three different heights, ground level, shade canopy, and sun canopy. 3. We implemented TreeTOP using two contrasting infrastructure concepts to demonstrate its applicability in both highly equipped canopy research facilities and forests without permanent canopy infrastructure. One installation relied on a canopy crane, grid power and fully automated irrigation, whereas the second was built by certified tree climbers and was equipped with an autonomous solar-powered, battery-operated irrigation system. At both sites, environmental sensor networks monitor the experiment. 4. TreeTOP successfully reproduced characteristic canopy microclimatic gradients, including increasing light availability, daytime air temperatures and thermal extremes with canopy height. Despite differing infrastructures, both implementations generated comparable microclimatic patterns, demonstrating that standardized canopy experiments are feasible in forests with or without permanent canopy access. By opening canopy space for manipulative experiments, TreeTOP provides a transferable framework for investigating plant performance, phenology, species interactions and microbiome assembly under realistic forest conditions.

Extrafloral nectaries (EFNs) are a widespread plant defense mutualism trait and are highly convergent, appearing in hundreds of plant lineages worldwide. Here we investigate a report of possible EFNs in Erythranthe angulosa, a recently-described California wildflower. We integrate field observations, insect bioassays, an induction experiment, and microscopy to test for signatures of EFN function, finding no evidence that the distinctive axillary swellings produced by E. angulosa function as EFNs. We also uncovered two distinct morphs at the type locality of E. angulosa that diverge in the number of axillary swellings produced, as well as other shoot architecture traits such as stem thickness, leaf size, and branch number. Although the axillary swellings appear to not function as EFNs, they remain a compelling morphological variant within the yellow monkeyflowers that may perform storage or another unknown function.

Urbanization is reshaping freshwater ecosystems, with well-documented effects across gradients of land-use change, hydrologic alteration, and habitat degradation. However, how biodiversity is organized among neighboring urban aquatic habitats that differ in hydrologic connectivity, disturbance transmission, residence time, management history, and opportunities for species movement is often less clear. This creates a challenge for interpreting urban fish communities at local scales as species occurrence may reflect both contemporary habitat filtering and historical contingencies including native persistence, interbasin transfer, stocking, and nonindigenous introductions. Here we use eDNA detections, historical records, phylogenetic information, and species trait data to investigate the fish assemblages of the Charlotte metropolitan region. We detect a highly mixed fauna that also depicts a strong signature of structured biodiversity profiles across taxonomic, phylogenetic, functional, and life-history dimensions between habitat types. In particular, bounded habitats contained assemblages with larger-bodied species that are fecund and faster to reproduce relative to free-flowing habitats. Species-level occurrence models did not support a simple trait-by-habitat rule. Instead our results demonstrate that urban aquatic habitats can sort historically mixed regional species pools into predictable assemblage-level life-history profiles while simultaneously retaining signatures of evolutionary and historical biogeographic contingency.

Surface phenotyping underpins plant science, preclinical animal research and entomology, yet across all three the measurement is almost always a photograph, which records a projection and not the surface itself. Here we present the Gentschinator3000, an open structured-light platform that brings high-end metric surface measurement within reach of laboratories with no optics expertise, combining documented open hardware, open reconstruction software and analysis workflows for under 4000 Euro in components. It resolves a planar reference to 45 m local flatness, registers full rotations to a loop closure of 156 m, and performs stably across acquisition ranges that we define. Applying one workflow to a leaf before and after desiccation, to murine anatomy and to a spread lepidopteran, we find that projection underestimates surface area by 11 to 41 %. That error grows with the condition under study, with the evaluation scale and with the direction of view, so it can confound phenotype comparisons dramatically. In murine limbs a 15-degree change of viewing direction shifts a projected inter-segment angle by up to 23.2 degrees, while the three-dimensional angle does not move. Projection geometry can therefore contribute as much to a measured phenotype as the biology it is meant to quantify.

A key goal in HIV-1 cure research is to understand why some individuals control viral rebound after stopping antiretroviral therapy (ART). Recent human studies have identified responding CD8+ T cells expressing Ki-67 and the transcription factor TCF-1 as correlates of post-treatment control, but the mechanistic basis of this association remains unclear. Using the theoretical framework of Conway and Perelson, we fit mechanistic within-host models to viral load and CD8+ T cell data from 9 individuals in a combination immunotherapy trial following ART interruption. Although Ki-67 and TCF-1 measurements were not used for fitting, the inferred effector cell expansion sensitivity, i.e., the responsiveness of effector expansion to low antigen levels, shows a strong linear relationship with Ki-67 and TCF-1 levels at rebound (Pearsons r {approx} 0.8). Building on this, we show analytically that the post-rebound viral load set point is inversely proportional to the effector cell expansion sensitivity, and thus strongly correlates with cycling (Ki-67+) CD8+ T cells (r {approx} -0.8) at rebound, and a subset that expresses TCF-1 (r {approx} -0.9). In effect, individuals with a larger proportion of CD8+ T cells responding to viral rebound, and a greater representation of TCF-1 expressing cells within the responding subset, achieve markedly lower viral set points through a higher effector cell expansion sensitivity. This mechanism is consistent with prior modeling in a non-intervention ATI setting, suggesting it may generalize across more rebound contexts. Our results provide a mechanistic explanation why both Ki-67+ responding CD8+ T cells and their TCF-1-expressing subset predict post-treatment control, linking clinical correlation to its underlying cause and highlighting Ki-67 and TCF-1 as potential early biomarkers of HIV immunotherapy success.

Most evidence of the consequences of climate change in natural systems has focussed on shifts in mean temperature (1,2), but the effects of extreme climatic events (ECEs) remain far less understood. This is particularly true for very severe ECEs that may occur only once every few decades. Understanding the consequences of these severe events for natural populations is nonetheless critical, since their frequency is predicted to rise under current climate change (3). Here we combine a unique long-term dataset spanning almost five decades of breeding (>20,000 events) and morphological data (>120,000 observations) in adult and nestling great tits (Parus major) and blue tits (Cyanistes caeruleus) with fine-scale temperature records to examine the effects of an unprecedented heatwave in May 2026 on breeding success and morphology. Average temperature during the heatwave (22-29 May 2026) was 7.85 C above the historical record, reaching +10.5 C (+4.32 SD) at its peak (25-26 May). These record-breaking temperatures significantly reduced adult breeding success and nestling bmass relative to expectation in the absence of a heat-wave. Given the heatwave was widespread (Fig. 1A), our findings from a single, exceptionally well-studied population are likely to generalise to other species exposed to the same event, providing key evidence that severe ECEs can substantially harm wild populations.

Standard methods for ancestral sequence reconstruction (ASR) rely on substitution models for the residues in a biological sequence and assume independent evolution across these sites, ignoring the epistatic interactions that shape molecular evolution. In contrast, deep learning models like variational autoencoders (VAEs) can learn low-dimensional representations ("embeddings") of sequences in a protein family that may implicitly handle these dependencies, raising the possibility of performing more accurate ASR by interpolating between extant sequence embeddings within the VAE's latent space. In this study, we test this hypothesis by developing and evaluating a VAE-based ASR pipeline. Benchmarking this approach against established likelihood-based and parsimony methods using various simulations of protein evolution, including scenarios with and without epistasis, we find that the VAE-based approach is consistently and significantly outperformed by standard methods, even in epistatic regimes where it was hypothesized to have an advantage. We further show that this failure is not due to a lack of phylogenetic structure in the latent space, which does contain evolutionary signal. Rather, the primary limitation is the information loss inherent to the autoencoding process: the VAE's decoder cannot generate sequences with sufficient fidelity for the precise demands of ASR.

A dental AI system and a dentist reading the same radiographs form a paired comparison. Published comparative studies often report the two arms separately against a reference standard, leaving the joint pattern of correctness between them unavailable for secondary paired inference. We show what that omission costs. The accuracy difference remains exactly identified; its sampling variance does not, so the report contains the estimate and not its uncertainty. On a study of 282 units, two published accuracies are consistent with 38 distinct joint tables whose confidence intervals differ in width by a factor of 2.5. The consequence is a three-zone decision map rather than a single threshold: differences at or below 1.06 points are non-significant under every compatible table, differences at or above 6.03 points are significant under every compatible table, and in between the published numbers cannot decide. We then show the omission is repairable at negligible cost. One additional integer, the number of units both arms classify correctly, identifies the joint table exactly and restores standard paired inference. For a panel of readers the pairwise dependences must arise from one joint distribution, a constraint that binds once three readers are present; publishing each reader's joint-correct count against a single reference reader cannot widen and may tighten every pairwise bound, and in a 7-arm experiment reduced them by a median of 37% even for pairs excluding that reference. Where the integer was never published we give DentalPair-Cert, an interval with finite-sample coverage uniformly over every admissible within-unit AI-dentist dependence under the independent-sampling-unit model, certified in both the nuisance maximization and the inversion. Across 4,200,000 simulated comparisons an independence analysis falls to 74.5% coverage with 12.2% type-I error; in a purposive sample of 9 recent comparative studies, 1 reported a paired test on discordant units.

Wildlife re-identification (re-ID) is a widely used and powerful tool with diverse applications in animal ecology and conservation. Current automated methods typically operate on single images of a single body part of the animal. However, a single encounter may contain multiple images capturing different body regions, each providing complementary individual-specific information. In contrast to automated approaches, researchers often manually select the most suitable images and regions for identification based on factors like visibility, occlusion and image quality. This creates a mismatch between automated methods and field practice, limiting the practical adoption of current automated re-ID pipelines. Here, we address this by introducing an encounter-based, multi-body-part re-ID framework, using sea turtles as a model taxon. Our framework combines three elements: (1) An orientation-aware deep learning model, TurtleDetector, that in addition to the full bodies, it also automatically segments key body regions, i.e. heads, front and hind flippers, from images within an encounter; (2) a hybrid body-part-specific retrieval method, that sequentially combines a fast global-feature model (MiewID or DINOv3) with a more accurate but costlier local-feature model (ALIKED with LightGlue); and (3) a merged identity-prediction strategy that selects the highest calibrated similarity score across all available body parts and images of an encounter. We evaluate the framework on three long-term re-ID datasets spanning three species, loggerheads, greens, and hawksbill turtles, under an evaluation protocol that mirrors real-world, time-aware re-ID workflows. Across datasets, combining multiple body regions consistently improved identification performance over the best-performing single body region, resulting to an increase of 4-6% in top-1 accuracy. Interestingly, body regions traditionally underused in sea turtle re-ID, such as the hind flippers and carapaces, provided complementary identifying information that improved encounter-level re-ID when integrated through the hybrid retrieval method. Our findings demonstrate that automated wildlife re-ID can benefit from moving beyond single-image, single-body-part identification towards encounter-level integration of all available visual evidence. Our work further suggests that, where feasible, field photo-acquisition protocols should aim to capture multiple informative views of an individual during each encounter. Importantly, many species and taxa, including elephants, primates, cetaceans, and other large vertebrates, possess such individual-specific features across multiple body regions, highlighting the broad potential applicability of our framework.

In vertebrates, top-down visual attention is a cognitive process where internal goals modulate the tuning of peripheral sensory systems. This leads to increased perceived contrast to both goal-relevant objects and areas of the visual field that are attended. Such a system would also be beneficial to bees, enabling them to detect and recognise the most profitable flowers in their environment. We tested whether bumblebees possess a top-down attentional system resembling that seen in vertebrates. We trained two groups of bees to collect rewards under high contrast targets. To potentially induce a difference in attention while searching for the targets, one group received a higher concentration of sucrose rewards compared to the other. During tests, the targets were presented with a series of lower contrasts to measure the contrast sensitivity curves of the bees induced by the different learnt reward levels. We predicted a stronger effect of any attention-like process on contrast sensitivity in the high reward group. We also repeated this experiment with the neonicotinoid pesticide imidacloprid dissolved in the sucrose rewards to test whether this affects bee attention. Across all test contrasts, higher rewards significantly increased bee accuracy when locating targets, lowered contrast thresholds and reduced the latency to make first choices. Imidacloprid reduced bee accuracy but did not influence first choice latency. These results suggest that learnt floral rewards can influence bee behavioural contrast sensitivity in a manner resembling vertebrate top-down attention and that imidacloprid may modulate this through effects on their nervous system.

Infections with bacteria which have evolved multi-drug resistance (MDR) cause millions of deaths worldwide. Large-scale efforts are gathering genotypic and phenotypic data on MDR bacteria, but methods for learning the structure, diversity, and predictors of evolutionary pathways to MDR have yet to take full advantage of these data. Here, we use evolutionary accumulation modelling (EvAM), an emerging class of machine learning methods with roots in cancer progression, to infer these evolutionary pathways across ESKAPEE pathogens (seven bacterial species that dominate health burdens), using a database of over 635k genotyped phenotypic observations from around the world. We identify global patterns in MDR evolutionary pathways, remarkably shared across multiple ESKAPEE species. Species-specific deviations from these stereotypical pathways are connected with geographical and demographic covariates, facilitating predictions of future MDR evolution. We verify these predictions with several hundred new phenotypes from ESKAPEE samples spanning decades of clinical infections in sub-Saharan Africa, demonstrating the capacity to forecast future MDR evolution from these inferred shared pathways.

Importance: Evidence on GLP-1/GIP therapies is largely derived from trials enrolling selected populations or medical records that miss utilization outside healthcare channels. No nationally representative cohort has characterized real-world uptake, indications, and access. Objective: To characterize GLP-1/GIP prevalence, indication, clinical profile, and access. Design: Prospective cohort study with three GLP-1/GIP surveillance waves (March 2024, December 2024, October 2025). Setting: The Understanding America Study, an address-based, nationally representative panel of approximately 15,000 US adults aged 18+ years initiated in 2014. Participants: UAS participants responding to at least one surveillance wave (n=9150). Exposures: GLP-1/GIP use status (never vs any use, comprising current and former use), self-reported primary indication (diabetes, weight loss, or other), and access pathway (traditional vs non-traditional). Main Outcomes and Measures: Survey-weighted prevalence of GLP-1/GIP use, overall and by indication and access pathway; sociodemographic, cardiometabolic, treatment, and access characteristics; and smartwatch-derived resting heart rate, heart rate variability, maximum activity heart rate, step count, and sleep duration and variability. Results: Among n=9150 adults (1274 with any use; 60.9% female; median age 53 years), weighted prevalence increased 46%, from 8.2% (March 2024) to 12.0% (October 2025) representing 32 million. Weight-loss indications grew, reaching nearly half of use (4.1% to 5.6%); diabetes-indicated use was stable (5.3% to 5.4%). Users carried high cardiometabolic burden (obesity, 68.2%; diabetes, 53.6%) but diverged by indication: diabetes-indicated users were older (median, 59 vs 49 years), whereas weight-loss-indicated users were more often female (69.9% vs 51.3%) and healthier. One in three users (~9 million) had non-traditional access, especially in weight-loss-indicated users, of whom 33% had no conventional prescription; 41% used compounding, online, or foreign pharmacies; and, 43% lacked coverage. Non-traditional users were five times as likely to report an unlisted, likely compounded formulation (19.8% vs 4.1%). All p<0.05. Conclusions and Relevance: Real-world GLP-1/GIP use has grown rapidly and diversified substantially in indication, access, and population profile. One in 3 users obtained treatment through nontraditional channels largely invisible to claims data, raising long-term safety, efficacy, and coverage questions. GLIMMER provides a public, nationally representative longitudinal evidence base for future payer and provider decisions.

Background & Aims: Fatigue is a central symptom of chronic liver disease (CLD), substantially impacting health-related quality of life (HRQoL). This study aimed to further understand CLD symptomatology, including fatigue, and its impact on HRQoL from a patient perspective. Methods: Abbott Global Assessment of Patients unmet needs (aGAP) was a multinational, cross-sectional survey in adults with compensated CLD in China, India and Mexico, conducted between July and November 2024. Adult participants who self-reported that they had physician-diagnosed CLD and were experiencing fatigue completed a quantitative survey to assess symptom burden and included three HRQoL patient-reported outcome (PRO) questionnaires (Patient-Reported Outcomes Measurement Information System [PROMIS]-29+2, Work Productivity and Activity Impairment - Specific Health Problem version 2.0 [WPAI: SHP], Multidimensional Fatigue Inventory [MFI]). Results: Overall, 505 participants (China: 200; Mexico: 105; India: 200) completed the study. Participants reported that their CLD-related fatigue sometimes, often or always affected their self-esteem/confidence (45.1%) and ability to maintain or acquire new employment (38.6%). Most participants reported moderate (51.3%) or serious (26.9%) fatigue, with 33.5% experiencing fatigue every day or almost every day. Many participants felt their social life was negatively impacted by their fatigue (47.3%) and that there were related financial difficulties (53.9%). Use of validated PRO tools demonstrated severe fatigue (MFI: overall mean [SD] 13.9 [3.4] general fatigue and 13.4 [3.6] physical fatigue) as well as substantial levels of work and activity impairment (WPAI: SHP overall mean [SD] 53.0 [26.4]) and high levels of anxiety, pain interference, depression and sleep interference (PROMIS T-scores [&ge;]54). Conclusions: Fatigue has a substantial impact on HRQoL among adults with CLD across several countries, highlighting a global unmet need for targeted interventions to effectively identify and manage the condition.

Muscle architecture shapes muscle function and changes with age, growth, training and disease, yet quantifying three-dimensional (3D) muscle architecture in vivo remains challenging. We introduce a hybrid fascicle tractography approach for freehand 3D ultrasound data that accurately reconstructs 3D muscle fascicles with respect to an objective, anatomically relevant coordinate system defined by the muscle's central aponeurosis. The hybrid approach combines Hessian-based fascicle detection with wavelet-based refinement to generate volumetric fascicle orientations. In a synthetic dataset with known ground truth, fascicle orientations and lengths were estimated with errors of [&le;]2{degrees} and ~1.5%, respectively. In vivo, the approach detected physiologically plausible fascicle lengthening in the human tibialis anterior following a passive plantar flexion rotation, whereas diffusion tensor imaging of the same muscle did not. The proposed method enables anatomically relevant, objective and non-invasive quantification of 3D muscle architecture in vivo, providing a practical framework for applications in clinical and applied muscle physiology.